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We numerically evaluate the entanglement spectrum (singular value decomposition of the wave- 
function) of paired states of fermions in two dimensions that break parity and time-reversal symme- 
tries, focusing on the spin-polarized px + ipy case. The entanglement spectrum of the weak-pairing 
(BCS) phase contains a Majorana zero mode, indicating non-Abelian topological order. In contrast, 
for the strong-pairing (BEC) phase, we find no such mode, consistent with Abelian topological order. 
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I. INTRODUCTION 

Two-dimensional fcrmion systems with pairing that 
breaks parity and time-reversal symmetries come in 
a variety of forms including quantum hall fluids,^ 
superfluids,^ superconductors,— and condensates of cold 
atoms near a Feshbach resonance^ For spin-polarized 
fermions, the simplest pairing order parameter that 
breaks these symmetries, Ap cc px + ipy, depends on 
the relative momentum p of the fermions in a pair. For 
momentum independent, s-wave pairing, a smooth cross- 
over occurs from weak-pairing (BCS) to strong-pairing 
(BEC). In the px + ipy case, the two phases have dif- 
ferent topological order and are separated by a quantum 
phase transitioni^ 

Recent proposals for fault-tolerant quantum computa- 
tion and information processing rely on topological order 
in fermion systems with p^ + ipy pairing,— but detect- 
ing and characterizing such order remain open problems. 
For example, the symmetry and bulk spectral proper- 
ties of the BCS and BEC phases are identical, but they 
have dramatically different topological order: quantum 
vortices have non-Abelian statistics in the weak-pairing 
phase and Abelian statistics in the strong-pairing phase. 
We apply ideas from quantum information to investigate 
topological order in these interesting, paired fermion sys- 
tems. 

The entanglement spectrum^ and the entanglement 
entropy^ contain information about the universal proper- 
ties of a quantum state. We define them by dividing the 
system into a block A with feature size L and an environ- 
ment B, and then performing a Schmidt decomposition, 

|V;)=5]e-^«-|V'.^)®|^f). (1) 

i 

Here, the orthonormal sets of states Ui'f)} span 

A and B. The entanglement spectrum {^i} gives the en- 
tanglement entropy S = J2i 

In this Letter, we report the first large-scale numerical 
calculations of the entanglement entropy and spectrum 
of two-dimensional fermion systems with px + ipy pair- 
ing. We find that the entanglement spectrum qualita- 



tively distinguishes the topological order occurring in the 
two phases. In particular, we find that the low-lying spec- 
trum in the weak-pairing phase contains a chiral, gapless 
fermion excitation. The weak-pairing phase is known to 
have a chiral, gapless Majorana edge modci^ This mode 
is related to the Majorana zero mode that appears in 
vortex cores and gives vortices non-Abelian statistics. 

We reduce the problem of evaluating the entanglement 
spectrum and entanglement entropy to diagonalizing a 
quadratic entanglement Hamiltonian.-'^^ This approach 
does not include fluctuations of the pairing order pa- 
rameter, and, hence, we do not expect to observe a uni- 
versal, topological term in the entanglement entropy^i 
in either the weak-pairing or strong-pairing phase?^ de- 
spite the fact that both phases have non-trivial quantum 
dimension D = 2. Indeed, we confirm that the size of 
the leading correction term depends on the geometry of 
the block, and is in fact proportional to the number of 
cornerSfiS, In contrast, the entanglement spectrum de- 
tects non-Abelian topological order in the ground-state 
wavefunction for states of paired fermions even when 
pairing fluctuations are neglected. 



II. PAIRING HAMILTONIAN 

The following BCS Hamiltoniar^ serves as a minimal 
model for a single band of spin-polarized fermions with 
Px + "i-Py pairing on a square lattice: 

H = (^-tclcr' - lr,r'clcl + h.C^ + 2A ^ C^C^ (2) 

{r,r') r 

We consider only nearest-neighbor (r, r') hopping t and 
pairing ^r,r' interactions. The hopping strength t and 
coupling A are taken to be real and positive, without 
loss of generality. The pairing interaction ^r,r' breaks 
both time-reversal and parity symmetries: 7r,r+x — 
—^r,r-x = i"fr,r+y = —ilr.i — y — il ■ Hcrc, 7 is real and 
X, y are the primitive translation vectors of the square 
lattice. We use periodic boundary conditions in our nu- 
merical calculations. 



FIG. 1: (Color online) Entanglement entropy S between a 
square of side length L and its environment as a function 
of A at fixed pairing strength 7 = 1.0. (Inset) The zero- 
temperature phase diagram of two-dimensional fermions with 
Px + ipy pairing and plot of the bulk spectral gap Eq. The 
phase boundary between weak-pairing and strong-pairing is 
the vertical 7-independent line at Ac = 2t. The spectral 
gap vanishes at the critical coupling and grows linearly with 
I A — Ac|. Data points indicate the parameters chosen in our 
numerical calculations (t=l). 

The pairing Hamiltonian ^ is quadratic and can 
be solved exactly using a Bogoliubov transformatioufi^ 
yielding the phase diagram shown in the inset of Fig. [T}^ 
The critical line at Ac = 2t, separates the weak-pairing 
(BCS) phase from the strong-pairing (BEC) phase. 
Both phases have a spectral gap Eq = t|A — Ad to 
bulk excitations show in the inset of Fig. [1] and de- 
termined by minimizing the Bogoliubov quasi-particle 

dispersion :i?p = \J + |Ap|^. The pairing order pa- 
rameter Ap = 27 (sin + isinpy) transforms under the 
symmetries of the square lattice in the same way as an 
i = 1,£^ = 1 spherical harmonic. At small p, we ex- 
pand Ap (X Px + ipy, and see the px + ipy pairing ex- 
plicitly. Similarly, at small p, the single-particle kinetic 
energy ^p = —2t{cospx + cospy) + 2A, takes the form 
^p ~ p^/2m* — /i, with effective mass m* = l/2t and 
fi — At — 2A. The weak-pairing phase A < Ac corresponds 
to /J, > 0, while strong-pairing A > Ac corresponds to 
/Lt < 0. Near the quantum phase transition fj, — 0, the 
low-energy spectrum Ep = ^ A.^'^p^ -V has a relativis- 
tic form with 27 playing the role of the speed of light. 



III. ENTANGLEMENT HAMILTONIAN 

The two-point correlation functions provide a complete 
description of the ground state of the quadratic Hamil- 
tonian (12) , and allow an efficient numerical evaluation of 
the Schmidt decomposition ([1]).^" In fact, the Schmidt de- 
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FIG. 2: (Color online) Leading correction term AS to the 
perimeter law for a square (a) and cross-shaped (b) parti- 
tion as a function of block size L. (Inset) geometry of the 
partitions. Notice the linear scale for AS and the logarith- 
mic scale for L in both (a) and (b). Solid lines are guides 
to the eye. Ratio AScr/ASsg of the leading correction terms 
from (a) and (b) as function of block size L: (c) within the 
weak-pairing and strong-pairing phases; (d) approaching the 
quantum phase transition from the strong-pairing regime. 

composition of the pairing Hamiltonian ground-state re- 
duces to diagonalizing the following entanglement Hamil- 
tonian _ffe which acts on the sites of the block A 

= ^ Cr,r' (4cr' + ^-c) + ^ {Frycic\, + /l.C.) .(3) 
r,r' r.r' 

Here, in contrast to the hopping parameters Cr.r' = 
/ (Pp/{2-Kfe^P<'-'''\Ep~£,p)/2Ep and pairing parame- 
ters F^,w = / d2p/(27r)2e'P (''"''')Ap/2£'p extend beyond 
nearest-neighbors and are given by the two-point corre- 
lation functions in the ground state of the pairing Hamil- 
tonian The entanglement Hamiltonian is quadratic, 
and can be exactly solved by numerically performing 
a Bogoliubov transformation to the quasi-particle op- 
erators a„, for n — ±1,±2,... ± Na, where, Na is 
the number of sites in the block A}^ In terms of the 
quasi-particles, the entanglement Hamiltonian has the 
form i/e = I]ti>o /(<^«)"nan> where, /(e) = (e" + 1)"^ 
is the Fermi function and the quasi-particle block en- 
ergies {e„} generate the entanglement spectrum. In 
particular, the entanglement entropy is given by 5 = 
-E„/(e")log/(e„). 

IV. RESULTS 

The entanglement entropy 5 as a function of the block 
size L is shown in Fig.[T] We consider various A sweeping 
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through the quantum phase transition, as shown in the 
inset. The entropy grows Unearly with L for this two- 
dimensional system. We interpret this as a perimeter law 
Sl = aL + . . . , where, the ratio of the correction terms to 
L vanishes in the limit L oo. Our large-scale numerical 
results agree with general arguments that a perimeter 
law must hold in the gapped phasesj-§ At the quantum 
critical point, the gap vanishes at a Majorana point,'' and 
no theoretical predictions or previous numerical results 
are available. 

Using these large-scale numerical results, we are able 
to extract the leading correction to the perimeter law 
AS* = —3(5' — aL)XL We plot the size dependence of 
the leading correction ASsq for the square shaped parti- 
tion shown in Fig. [5)^ a) and for the cross-shaped partition 
AS'cr shown in Fig.O^b). For both geometries, the lead- 
ing correction grows at the critical point with L, without 
sign of saturation. By contrast, in the weak-pairing and 
strong-pairing phases, the leading correction saturates to 
an L independent value as L oo. We interpret the 
growth at the critical point as a logarithmic divergence, 
of the form S — aL — b log L + . . . . This is the first indica- 
tion that a Majorana point exhibits a logarithmic correc- 
tion to the perimeter law for the entanglement entropy, 
although additive logarithmic corrections have been ob- 
served in two-dimensional systems with other kinds of 
nodal excitationSfiiii^ 

In Fig. ^c) and (d), we analyze the geometry de- 
pendence by plotting the ratio of the leading correction 
AScr/ ^Ssq for the two partition geometries. In both 
strong-pairing and weak-pairing phases (Fig. (2};), the ra- 
tio AScr/ ASsq — > 3 approaches the ratio of the number 
of corners in the cross partition to the number in the 
square partition. We have examined other geometries 
and find the behavior AS = cUc, where, ric is the num- 
ber of corners and c is a positive coefficient.— In con- 
trast, when pairing fluctuations are allowed, the topolog- 
ical term AS' = 3 log 2 has no geometry dependence.^^ 
To check that our results reflect the asymptotic behavior 
of the system, we reduce the de-tuning from the criti- 
cal point while staying on the BEG side (See Fig. 2(d)). 
For system size L » ^ much bigger than the diverg- 
ing length scale ^ — 2'^/\fi\ oc |A — Ac|~^, the behavior 
AScr/ ASsq — > 3 observed deep within the gapped phases 
(Fig. 2(c)) emerges near the critical point as well. Thus, 
these large-scale numerical simulations indicate a geo- 
metric origin of the leading correction to the perimeter 
law for the entanglement entropy. 

To detect topological order, we turn to the entangle- 
ment spectrum shown in Fig. [31 Now, in the weak- 
pairing phase, the energy spectrum of the pairing Hamil- 
tonian ^ for a system in the form of a disc of ra- 
dius R contains a chiral fcrmion edge mode with energy 
E cx m/R proportional to angular momentum To 
detect such a mode in the square geometry, one must la- 
bel the quasi-particle block energies {e„} by the phase 
factor (/)„ =0,7r/2,7r,37r/2 acquired by the quasi-particle 
wavefunction under the elementary n/2 rotation symme- 
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FIG. 3: Low-lying quasi-particle entanglement spectrum 
(a) in the weak-pairing phase, (b) at the quantum phase tran- 
sition, and (c) in the strong-pairing phase with fixed pairing 
strength 7 = 1.0 and system size L = 24. We divide the spec- 
trum into four sectors, corresponding to the irreducible repre- 
sentations of the point group of the square lattice and labeled 
by the phase factor acquired by the quasi-particle wavefunc- 
tion during a 7r/2 rotation. 



try of the square lattice. This phase factor plays the role 
of angular momentum in a lattice system. 

In the weak-pairing phase. Fig. [3fa) we find that both 
the energy e„ oc n and the phase factor 20„/7r = n(mod4) 
are proportional to the level index n — 1,2,.... Elimi- 
nating the level index, we find en oc (pn- This is precisely 
the relationship expected for a gapless chiral mode, and 
observed in the weak-pairing phase of the pairing Hamil- 
tonian ^ along an edgci^ In the strong-pairing phase. 
Fig. [SKc), the phase factor 0„ and level index have no 
apparent relationship. At the critical point Fig. ^h) 
the phase factor and level index are proportional for the 
lowest levels, but have no relationship for higher quasi- 
particle block energy. We contrast the dispersionless low- 
lying spectrum in the strong-pairing phase with the lin- 
early dispersing spectrum in the weak-pairing phase, and 
compare the weak-pairing result e oc to the energy 
spectrum i? oc to of the pairing Hamiltonian ([2]) in the 
weak-pairing phase. 

To test the identification further, we show in Fig. 21 
the finite-size scaling of the minimum quasi-particle block 
energy ei plotted on a log-log scale at fixed pairing ampli- 
tude 7 = 1.0. For the pairing Hamiltonian (2) on a disc 
of radius R, the minimum quasi-particle energy scales as 
El (X 1 /R in the weak-pairing phase and tends to a con- 
stant in the strong-pairing phase, as i? — > 00. In Fig. 21 
the data for the strong-pairing phase A > Ac tend to a 
constant as L — )■ 00. By contrast, in the weak-pairing 
phase A < 2.0, the minimum block energy drops to zero 
ei ~ 1/L, for system sizes L 3> ^ large compared to 
the diverging length scale ^ = 27/IA — Ad characteriz- 
ing critical fiuctuations. In the quantum critical regime. 
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FIG. 4: (Color online) Finite-size scaling of the minimum 
quasi-particle block energy ei plotted on a log-log scale at 
fixed pairing amplitude 7 = f .0. In the weak-pairing phase 
A < 2.0, the dashed lines are best fits to the scaling form 
ei ~ 1/L. 



L <^ ^, the finite-size scaling of the minimum quasi- 
particle energy is intermediate between those of weak 
and strong-pairing phases. Remarkably, the contrast in 
Fig. 4 between the finite-size scaling of the weak-pairing 



and strong-pairing phases occurs even for relatively small 
block sizes L/£^ « 1. On the other hand, the data in Fig. 
2(c) and (d) show that the finite-size corrections to the 
entanglement entropy require significantly larger systems 
L/^ « 3 to see the asymptotic behavior. 



V. CONCLUSION 

In this Letter, we study topological order in paired 
states of fermions with parity and time-reversal symme- 
try breaking. Large-scale numerical calculations of the 
entanglement spectrum and entanglement entropy reveal 
universal behavior. In particular, we find a chiral, gapless 
Majorana fermion excitation in the entanglement spec- 
trum of the weak-pairing phase, and contrast this with 
the gapped spectrum in the strong-pairing phase. A va- 
riety of topological phases can be described by a pairing 
Hamiltonian that neglects order parameter fluctuations. 
We suggest that large-scale numerical calculations of the 
entanglement spectrum are a robust way to detect and 
characterize non-Abelian topological order in the ground- 
state wavefunction of such phases. 
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